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Abstract. We investigate gas dynamics in the presence of a double inner Lindblad 
resonance within a barred disk galaxy. Using an example of a prominent spiral, MlOO, 
we reproduce the basic central morphology, including four dominant regions of star 
formation corresponding to the compression maxima in the gas. These active star 
forming sites delineate an inner boundary (so-called nuclear ring) of a rather broad 
oval detected in the near-infrared. We find that inclusion of self-gravitational effects in 
the gas is necessary in order to understand its behavior in the vicinity of the resonances 
and its subsequent evolution. The self-gravity of the gas is also crucial to estimate the 
effect of a massive nuclear ring on periodic orbits in the stellar bar. 
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^ Invited talk at the Centennial Nobel Symposium on Barred Galaxies and Circum- 
nuclear Activity, A. Sandquist et al. (Eds.), Springer- Vcrlag, 1996 
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1 Nuclear Starburst Galaxies 

The role of active galaxies within the framework of galactic evolution is far from 
clear. In particular, it is still unknown if central activity encompasses a small 
percentage of galaxies or if it is a normal evolutionary stage. It also remains 
undetermined whether two notable types of such activity, nuclear starbursts and 
Seyfert nuclei, are 'genetically' linked. 

A subgroup of disk galaxies shows intense star forming activity within the 
central few hundred parsccs, in the so-called nuclear rings. In the visual and 
ultraviolet wavelength ranges and at a high spatial resolution, these rings fre- 
quently appear to be patchy and incomplete, and/or consist of a pair of tightly 
wound spiral arms (Buta and Crocker 1993). In the near- infrared (NIR), they are 
regular and weakly elliptical (Knapen et al. 1995a; Shaw et al. 1995). Molecular 
gas distribution based on CO emission reveals a complex morphology at and in- 
terior to the rings (review by Kenney 1996). The origin of these rings is related 
to the inner Lindblad resonance (ILR), i.e. the resonance between the planar 
stellar orbits and the perturbing force of a stellar bar or of an oval distortion. 
Although self-gravity in the gas was ignored in earlier numerical simulations of 
disk galaxies, these simulations clearly showed that gas is focused into nearly 
circular orbits interior to the outer ILR (OILR; e.g. Schwarz 1984; Combes and 
Gerin 1985). As the gas accumulates in this region, no further evolution occurs 
and it was suggested that the gas is converted into stars with a high efficiency 
(Elmegreen 1994). More sophisticated 2D modeling involving self-gravity in the 
gas supported this picture (Shaw et al. 1993), indicating at the same time that 
increasing surface density in the gas may lead to its fragmentation and destruc- 
tion of the gaseous ring (Wada and Habe 1992). Effects of fragmentation are 
lessened if energy deposition by massive stars is taken into account by inducing 
turbulence in the gas (Heller and Shlosman 1994). 

Recent observations of molecular gas distribution in the centers of a num- 
ber of barred galaxies provide some indirect evidence (based on stellar rotation 
curves) that a double ILR may in fact reside there (Kenney et al. 1992; Knapen 
et al. 1995a). Although gas flows in barred galaxies have been modeled for about 
two decades now, gas behavior in the vicinity of a double ILR is not well under- 
stood. Due to the low resolution of numerical schemes and to a common belief 
that the inner ILR (IILR) is probably located too close to the center to be obser- 
vationally resolved, most efforts have been devoted to the study of gas dynamics 
between the OILR and outer Lindblad resonance (OLR; e.g. Athanassoula 1992). 
For these reasons, numerical schemes usually fail to catch the nuclear ring phe- 
nomenon. 

Here we investigate a self-consistent gas evolution in the central resonance 
region of a moderately barred galaxy. Effects of star formation are incorporated 
at some level. For convenience, we choose the MlOO-like total mass distribution 
which was claimed to possess a double ILR (Arsenault et al. 1988; Knapen et 
al. 1995a). Furthermore, we highlight the effect that massive nuclear rings have 
on the dominant stellar orbits in the galactic disk, i.e. the back-reaction of the 
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stellar component to the self-gravitating gas accumulating between the ILRs. A 
full account can be found in Knapen et al. (1995b) and Heller and Shlosman 
(1996). 

2 Twisting of NIR Isophotes in MlOO 

MlOO, the brightest (barred) spiral galaxy in Virgo, displays all the virtues of 
a nuclear starburst and is inclined at 30° ± 3°. Surprisingly, the UV/optical 
starburst ring formed by a tight pair of spirals, lies around 10" from the center 
(1" corresponds to 83 pc at a distance of 17.1 Mpc), whereas its NIR (2.2/tm) 
counterpart has a substantially larger radius of ^ 18" and a considerable width 
of ~ 16" (Knapen et al. 1995a). It is weakly oval (minimal elhpticity ~ 0.13) 
and its semimajor axis leads the stehar bar by about 60° — 70°. Most unusual is 
the observed gradual twist of the NIR isophotes, from the position angle (P.A.) 
of the ring towards the P.A. of the stellar bar, both exterior and interior to 
the ring (see Fig. 1). In other words, the large-scale bar is dissected by an oval 
ringlike zone oblique to the bar. The existence of this inner bar-like feature 
was confirmed independently by high-resolution CO observations of the velocity 
field (Rand et al. 1996). Molecular gas within the central 10" participates in high 
velocity ~ 100 km s^^ non-circular motions along the P.A. of the stellar bar. 



Fig.l. The K and I iso- 
photes in MlOO: position 
angle (top) and elliptici- 

ty {bottom) (From Kna- 
pen et al. 1995b). 



Stehar orbits are expected to change their orientation abruptly by 90° at 
each resonance in the disk, which can be understood within the framework of 
forced oscillations. These orbits arc oriented along the stellar bar between the 
OILR and the corotation radius (xi family), they are perpendicular to the bar 
between the IILR and OILR {x2 family), and are again aligned with the bar 
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between the center and the IILR. NIR isophotes in a quiescent galaxy are be- 
Heved to follow the overall mass distribution, i.e. they arise from an old stellar 
population represented by K and M giants, (e.g. Frogel 1985). However, gradual 
skewing of isophotes in MlOO is indicative of gas behavior across the resonances 
rather than the behavior of stars. There are two possible ways to explain this 
twisting of NIR isophotes between 10" to 30" in MlOO. First, there can be a 
sufficient contribution from massive young stars, K and M supcrgiants (Knapcn 
et al. 1995a). These stars arc expected to be found in the star bursting regions, 
and, being dynamically young (age less than 10''^ yrs), they should follow gas 
rather than stellar orbits. In such a case, the NIR light will not follow the mass 
distribution. Knapen et al. (1995b) found that this is the most plausible explana- 
tion for two NIR 'hot spots' in MlOO, and similar conclusions have been reached 
for active star forming regions in NGC 1309 (Rhoads 1996). Alternatively, gas 
gravity can drag some of the old population stars in the stellar bar towards the 
X2 orbits in the ring (Shaw et al. 1993). As we show in Section 4, the growing 
accumulation of gas between the ILRs is capable of affecting the main periodic 
orbits in the bar (especially the Xi orbits outside the OILR) in such a way as 
to support a gradual twist of stellar orbits, from being aligned with the bar to 
becoming almost perpendicular to it. 



3 Simulations of Stellar and Gas Dynamics in MlOO 

To further understand the circumnuclear morphology in MlOO and to confirm 
that it is compatible with the presense of a double ILR there, we have tailored our 
numerical simulations to a M100-/«A;e mass distribution. Stars and gas, embedded 
in halo and bulge potentials, were evolved by means of a 3D hybrid SFH/N- 
body code (Heller and Shlosman 1994). Models without and with star formation 
have been constructed, Ql and Q2 respectively. In the Q2 model, the gas is 
considered to undergo 'star formation' if it is Jeans unstable, if it participates 
in a locally converging flow, and if its density exceeds 20 Mq pc^"'. We assume 
that only massive OB stars form, which deposit 10^^ ergs per 10^ yrs per star 
in the gas by means of their line-driven winds, as well as deposit 10^^ ergs 
per 10^ yrs as supernovae, leaving no remnants. This energy is assumed to be 
instantly radiated away by the gas, with only 5% being retained and converted 
into turbulent motion. One unit of time, [r], corresponds to 3.75 x 10^ yrs. 

Initially, an axisymmetric model of a galactic disk was constructed which was 
dynamically unstable and developed a stellar bar in the process of evolution. The 
resulting bar strength was moderate, q ~ 0.3, where q is deflned as a maximum 
ratio of the sum of m = 2,4,6, and 8 Fourier components of the nonaxisym- 
metric gravitational force to the m = axisymmetric component. Based on the 
(axisymmetric) rotation curve of the model after ~ 3 rotation periods, a double 
ILR could be found in the central region, IILR at ~ 240 pc and OILR at ~ 1.4 
kpc. However, the epicyclic approximation is invalid because of the relatively 
strong bar. Instead, we have analyzed the dominant families of periodic orbits 
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in the model potential. It is customary to extend the ILRs into the nonlinear 
regime based on the radial extent of the X2 family of orbits. The corresponding 
limiting semimajor orbital axes in Ql and Q2 have been found at ~ 500 pc and 
~ 1.3 kpc, and adopted as the nonlinear IILR and OILR, respectively. Hence, 
the resonance region between the ILRs was reduced substantially compared to 
the linear regime. 

After an initial transient in the Ql model, the gas formed a pair of trailing 
shocks along the leading edge of the stellar bar, in agreement with 2D numerical 
simulations (e.g. Athanassoula 1992). Deeper in the potential well, gas dynam- 
ics in the vicinity of ILRs was dominated by a pair of tightly wound trailing 
shocks and a pair of leading shocks (see Fig. 2, at r = 15). Such a shock system 
in the vicinity of a double ILR is predicted on the base of a simple epicyclic 
approximation which is not adequate at this stage. In particular, we note that 
two systems of shocks, trailing and leading, interact non-linearly through a pro- 
nounced cuspy feature (caustic; Fig. 2, r ^ 17). To the extent that this shock 
system delineates spiral arms, we observe a transient pseudo-ring made out of 
a pair of tightly wound spirals between the ILRs. In fact, two separate gas cir- 
culations (oval 'rings') form between the ILRs due to the action of the shocks. 
The outer gaseous ring evolves as to (almost) align itself with the minor axis of 
the stellar bar, while the inner gaseous ring (almost) aligns itself with the major 
axis of the bar. 

This evolution comes about as a response to the gravitational torques by the 
bar and in order to minimize them. However, the effects of self-gravity in the gas 
are crucial in order to understand its dynamics. The self- gravity acts in a manner 
similar to surface tension: the outer gaseous ring settles down on lower energy 
X2 orbits, away from the OILR. At the same time, the inner ring shrinks across 
the IILR and settles on the xi orbits. This behavior is depicted in Figure 2, at 
around r ~ 17 — 24. It is imperative to mention that the gas is not locked between 
the ILRs, but merely experiences a temporary slowdown there. The outer ring is 
constantly perturbed by a number of density inhomogeneities and 'rains' onto the 
inner ring which contains orbits deeper in the potential well. (In the presence 
of star formation, as we discuss below, these perturbations are caused by the 
turbulent motions in the gas excited by stellar winds and supernovae.) As the 
inner ring becomes more massive, it is subject to gravitational fragmentation. 
Neglecting the effects of star formation at this stage is not justified. 

Technically speaking, star formation in the Q2 model has introduced 'tur- 
bulent' motion in the gas and induced mixing between material with a different 
angular momentum and energy. As a result, gaseous circulations between the 
ILRs widen notably and merge, and it is more appropriate to speak of a gaseous 
disk which extends inwards, from the outermost X2 orbits, much of the way to 
the center. This disk is an oval, and its semimajor axis is positioned almost at 
right angles to the bar. 

As expected, the shock system outlines the regions of intense star formation 
in the Q2 model (see Fig. 3). However, there is no one-to-one correspondence 
between the shocked and gravitationally unstable gas, as the shock and star 



6 Isaac Shlosman 



Fig. 2. Logarithmic gray-scale map of shock dissipation inside the OILR in MlOO (Ql 
model, without star formation). Time is given in the upper left corners, [r = 1] = 
3.75 X lO'^ yrs. Each frame is 2.6 kpc across. The gas flows counterclockwise, the stellar 
bar is horizontal. (Prom Knapen et al. 1995b). 

formation maps show. At the same time, four major sites of star formation in the 
resonance region persist during most of the simulation time — all corresponding 
to the maxima of dissipation in the gas. Two elongated star forming regions 
are found downstream from the place where the inflow along the outer shocks 
crosses the bar's minor axis and encounters the gas circulation on the X2 orbits. 
Azimuthal smearing of star forming regions is a direct consequence of the time 
scale for Jeans instability becoming an appreciable fraction of the orbital time 
scale so close to the rotation axis. Another pair of star forming regions is located 
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around the IILR (r ~ 10 — 25, in Fig. 3), slightly ahead of the bar's major axis 
and where the cuspy feature is seen in the Ql model. This prevailing morphology 
dominated by four star forming regions (at and just outside the IILR) appears 
to be robust during most of the simulation time. It is compatible with the loci of 
star formation in the U, V, and Ha images of the inner 2 kpc in MlOO (Knapen 
et al. 1995a). No star formation correlates with the position of the OILR. 

The star formation rate in the Q2 model reaches its peak around r ~ 25 — 28, 
exhibiting burst behavior with a typical time scale of ~ 10^ yrs. Around r ~ 28, 
the mass inflow rate across the IILR peaks strongly, indicating a catastrophic loss 
of angular momentum by the gas (r = 28.4, Fig. 3). The subsequent evolution of 
this dynamical runaway in the self-gravitating gas was discussed by Heller and 
Shlosman (1994) in the absence of star formation. We view this process as the 
initial phase of decoupling of the gaseous bar from the large-scale stellar bar, 
as envisioned in the 'bars within bars' scenario (Shlosman, Frank and Begelman 
1989). It is not clear how much of the underlying old population participates in 
this process (and if it does at all) . The amount of gas at the onset of instability, 
~fewx 10® Mq, is ~ 10% — 20% of the mass interior to the runaway region, 
which is more than an order of magnitude less than quoted by Ho, Filippenko 
and Sargent (1996). This amount of gas is available even in the early-type disks. 

In the more realistic case, advanced here in the Q2 model, the outcome of 
evolution depends in a sensitive way on the efSciency of star formation, the 
state of the interstellar medium, and its ability to retain energy deposited by 
massive stars. These questions are potentially relevant for our understanding of 
fueling the nonthermal activity in Seyfert nuclei and must be addressed in future 
work. We conclude, that a characteristic time for the gas to 'filter' through the 
resonance region is ^ 10^ yrs. This can be taken as a rough estimate of a nuclear 
ring's lifetime. 

4 Massive Nuclear Rings: Affecting Periodic Orbits 

The above numerical simulation of gas dynamics in the presence of a double 
ILR demonstrates explicitly that the radial gas inflow is not stopped between 
the resonances but only slows down there — a kind of 'self-mulching lawn mower' 
effect. This results in gas accumulation in the form of a massive elliptical 'ring'. 
Molecular mass of as much as a fewxlO^MQ stored in the ring is not out of 
the question. Besides enhanced star formation caused by favorable conditions in 
the region, such a ring will have gravitational effects on the stellar component 
in the disk. Here, we are mainly interested in how the main periodic orbits are 
affected, and what consequences it may have on the gas circulation in the bar. 

For simplicity we use a 3D analytical model of a galaxy consisting of a disk, 
bulge, halo and stellar bar. Furthermore, in an attempt to simulate the effect 
of a massive nuclear ring within the central kpc, we make assumptions about 
its shape and the orientation of its major axis with respect to the stellar bar. 
Results of a 3D orbit analysis for such gravitational potentials are presented 
elsewhere (Heller and Shlosman 1996). Here we discuss only necessary details. 
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Fig. 3. Left: Logarithmic gray-scale map of shock dissipation inside the OILR in MlOO 
(model with star formation). Time is given in the upper left corners, [t = 1] = 3.75 x lO'^ 
yrs. Each frame is 2.6 kpc across. The gas flows counterclockwise, the stellar bar is 
horizontal. Right: Star formation map corresponding to the region shown on the left. 
(Prom Knapen et al. 1995b). 



The most pronounced change in the periodic orbits, when a circular ring 
(or elhptical ring whose major axis coincides with that of the bar) is added, is 
that the extent of the ILR resonance region is increased with the ring's mass, 
weakening the stellar bar. In addition, the xi orbits of different energy intersect 
in the vicinity of the ring (as do X2 orbits within the ring) . This has a two-fold 
effect on the gas: the phase space available to X2 has increased, but at the same 
time, orbits near the ring became intersecting and unable to hold gas, amplifying 
dissipation there. Thus, we expect that the gas will fall through the IILR after an 
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initial stage of accumulation, exactly as observed in the numerical simulations. 

An additional and qualitatively different effect is obtained when the ring is 
mildly elliptical and its major axis is oblique to the bar, leading it in the direction 
of galactic rotation. This configuration is the one typically observed in nuclear 
starburst galaxies and appears as a long-lived transient in our numerical simula- 
tions. The behavior of xi orbits can be qualitatively understood in this case as a 
response to the perturbing forces of the stellar bar and of the oblique ring. Both 
forces have the same driving frequency but are phase-shifted. A straightforward 
application of an epicyclic approximation to the motion of a viscous 'fluid' re- 
veals the rich variety of possible responses in the gas to this driving force. A 
representative case, calculated using fully nonlinear orbit analysis when the ring 
leads the stellar bar by 60°, is shown in Figure 4. The change in the position 
angle of the Xi with distance to the ring is rather dramatic, starting at large 
radii with a slowly growing phase shift which reaches a maximum of 34° in the 
leading direction, followed by a rapid decline to 0° just interior to the ring and 
then continuing to -11° in the trailing direction deep into the bulge. 



Fig. 4. Twisting of xi periodic orbits supporting the stellar bar in the presence of 
a massive 10^ Mq oblique elliptical ring (e = 0.4). Ring's potential is softened with 
e = 100 pc. The stellar bar is horizontal. Rotation is counterclockwise and corotation 
is at 5 kpc. The ring's semimajor axis ([r = 0.04] = 400 pc) is offset to the bar by 
60° in the leading direction. The frame is 2 kpc on a side. (From Heller and Shlosman 
1996). 



5 Implications: Dynamics of the Circumnuclear Region 



Our dynamical study of gas flow across a double ILR resonance in the circum- 
nuclear region of a barred galaxy provides some insight into relevant processes 
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which accompany this flow. First, self-gravitating effects in the gas accumulat- 
ing between the resonances are crucial in understanding its subsequent evolution 
(besides the star formation). In particular, self-gravity acts as a 'surface tension' 
and the gas moves deeper into the potential well, away from the OILR, and across 
the IILR. Star formation is clearly peaked around the IILR, partly because the 
gravitational torque's sign changes across this resonance causing additional com- 
pression in the gas and creating conditions favorable to Jeans instability there. 

Second, we flnd a complicated but basically a ring-like morphology for the 
distribution of star forming regions in the nuclear starbursts. Dominant regions 
of star formation lie downstream from two main compression sites of the gas 
('twin peaks'), when it crosses the minor axis of the stellar bar. This is somewhat 
outside the IILR. Additional pair of star forming regions are found at the IILR, 
on the major axis of the bar and slightly leading it, and correspond to the twin 
'hot spots' found in MIOO. 

Third, there is no prominent star formation associated with the OILR. This 
explains why the UV/optical and Ha observations identify the star forming ring 
in MIOO inside its NIR ring. Unless the star formation, already at this stage, is 
very efficient in consuming the gas at the IILR, the gas 'breaks through' and 
falls towards the center, an event that is accompanied by a prominent burst of 
star formation. This dynamical runaway of self-gravitating gas inside the IILR 
depends on details of star formation and physics of the interstellar medium, and 
it is outside the scope of this study. We only comment that the characteristic 
time scale of this runaway is ~fewxlO'' yrs, much shorter than the time it takes 
for the gas to 'filter' through the resonance region, ^ 10^ yrs. 

Fourth, we find that massive nuclear rings are capable of perturbing the 
gravitational potential in the circumnuclear regions, thus affecting the main 
periodic orbits there. The phase space allowed to the orbits aligned with the 
minor axis of the bar {x2 family) is substantially increased and the orbits aligned 
with the bar (xi family) are significantly distorted. Orbits with different values of 
the Jacobi integral are found to intersect, meaning that gas cannot be retained 
there and will move inwards across the IILR. In the most interesting case of 
an elliptical ring oblique to the bar, Xi orbits are gradually twisted, in a way 
similar to the skewing of NIR isophotes observed in MIOO. So these orbits can 
trap both gas (accompanied by K and M supergiants) and old population stars. 
It may be possible, in a such a case, to estimate the 'degree of asymmetry' in 
the nuclear ring (i.e. azimuthal mass distribution) based on the observed change 
in the ellipticity and position angle of NIR isophotes with radius. 

To summarize, the gas seems to be the prime dynamic agent in the circum- 
nuclear regions of at least some disk galaxies, although its mass is only a fraction 
of the mass in the stellar component there. Further observations of molecular 
gas kinematics in active and normal galaxies will provide clues to understanding 
their central activity. 
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